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L SUMMARY 

»/l  submariner  is  able  to  estimate  the  course  of  a ship  from  his  record  of  successive 
bearings  on  the  ship.  The  “ bearings  only ” technique  provides  a mathematical  approximation 
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INTRODUCTION 


Passive  listening  devices  available  to  the  submariner  enable  him  to  construct  a record  of 
the  successive  bearings  of  a ship  or  convoy,  though  they  may  provide  little  or  no  information 
on  the  range  of  the  potential  target.  A mathematical  approximation  known  as  the  “bearings 
only”  approximation  provides  an  estimate  for  the  angle  to  bow  (ATB)  of  the  convoy.  The 
mathematical  formula  used  does  not  require  any  information  about  target  range,  nor  does  it 
provide  any  such  information.  However,  a knowledge  of  the  maximum  capabilities  or  usual 
operating  procedures  on  potential  targets,  or  the  intensity  of  sound  received,  may  enable  him 
to  guess  an  approximate  range. 

An  estimated  ATB  obtained  from  the  “bearings  only”  formula  will  be  subject  to  three 
sources  of  error: 

(1)  Non-uniform  motion  of  the  target 

(2)  Errors  in  the  measurement  of  bearings 

(3)  Error  introduced  by  the  approximation  involved  in  the  calculations. 


This  report  records  the  results  of  an  investigation  of  the  third  type  of  error  above.  A large 
number  of  representative  situations  involving  a submarine  and  a convoy  have  been  examined, 
and  the  estimated  ATB  compared  with  the  known  ATB  in  each  case.  The  data  generated  here 
provide  an  indication  as  to  when  (in  the  absence  of  the  first  two  types  of  error  listed  above) 
the  estimated  ATB  from  the  “bearings  only”  approach  is  most  likely  to  be  reliable. 


Time  Bearing  Plot 


The  submariner  plots  the  received  bearings  of  a target  as  a function  of  time  on  a “time 
bearing  plot”.  A smooth  curve  is  drawn  through  these  plotted  points.  All  of  the  submariner’s 
calculations  involve  bearings  read  off  this  smooth  curve  rather  than  the  original  raw  data  plotted. 
This  is  done  in  an  attempt  to  minimize  the  effect  of  bearing  measurement  errors. 

In  addition,  if  sufficient  time  is  available,  the  submariner  may  develop  time  bearing  plots 
for  two  submarine  courses;  one  directed  along,  and  the  other  approximately  at  right  angles  to, 
the  radius  vector  to  the  target. 

Estimation  of  Convoy  Course 

In  making  his  estimate  of  target  motion,  the  submariner  assumes  that  the  convoy  maintains 
constant  speed  and  direction  while  he  is  taking  bearings  on  it.  From  the  time  bearing  curve 
the  submariner  interpolates  a series  of  bearings  corresponding  to  equal  time  intervals  between 
adjacent  bearings.  These  bearings  are  plotted,  and  a line  is  drawn  across  them  to  give  equal 
intercepts  travelled  by  the  convoy  in  equal  time  intervals.  The  direction  of  this  line  then  gives 
the  approximate  course  of  the  convoy. 

Suppose  that  three  such  bearing  lines  are  drawn;  bearing  B at  the  end  of  the  series,  bearing 
C in  the  middle  and  bearing  E at  the  beginning.  Allowing  for  the  motion  of  the  submarine, 
there  are  three  different  regions  in  which  the  convoy’s  course  might  lie  for  any  given  B , C and  E. 
Three  possible  courses,  one  in  each  region,  are  shown  in  figure  1.  Normally  the  submarine’s 
detection  equipment  will  be  able  to  distinguish  between  course  I (figure  1)  and  the  other  two. 
Additional  information  would  be  required  to  distinguish  between  courses  2 and  3. 

The  geometrical  situation  corresponding  to  convoy  course  I is  illustrated  in  figure  2. 
Similarly,  course  2 is  illustrated  in  figure  4.  In  figures  2 and  4 the  ATB  at  the  time  corresponding 
to  bearing  B is  A o.  During  the  time  interval  between  bearings  E and  B the  convoy  moves  from 
(x2,  y2)  through  (.v„  v,)  to  (O,  P)  while  the  submarine  moves  from  (x4,  yt)  through  (x3,  y3)  to 
(0,  0).  This  geometrical  layout  corresponds  to  the  location  of  cartesian  axes  with  the  origin  at 
the  final  submarine  position,  and  the  >’-axis  passing  through  the  final  position  of  the  convoy. 
Any  other  cartesian  co-ordinate  system  may  be  made  consistent  with  this  scheme  by  suitable 
translation  and/or  rotation  of  co-ordinate  axes. 
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The  “bearings  only"  approximation  for  ATB  rests  on  two  main  assumptions: 

(1)  The  convoy  maintains  a straight  course  and  moves  equal  distances  in  equal  time  intervals. 

(2)  The  motion  of  the  submarine  may  be  neglected. 

All  of  the  possible  relative  orientations  of  submarine  and  convoy  are  represented  by  figures  2 
and  4.  Application  of  the  above  assumptions  to  figures  2 and  4 allows  their  simplification  to 
figures  3 and  5 respectively.  Each  case  then  yields  a unique  solution,  Ev  for  the  approximate 
ATB.  This  solution  is  derived  in  Appendix  1. 

THE  CALCULATIONS 

In  figures  2 and  4 the  submarine  moves  a distance  of  W units  in  one  unit  of  time.  In  the 
same  time  the  convoy  moves  V units.  The  submarine  and  convoy  are  P units  apart  when  the 
final  bearing  is  recorded.  Substituting  W — 1 in  the  calculation  allows  P and  I'  to  be  used  as 
ratios,  thus  obtaining  maximum  generality  in  the  results. 

The  ratios  chosen  for  the  calculations  were  based  on  the  following  data: 

Submarine  speed  = 5 to  15  knots 
Convoy  speed  = 10  to  25  knots 
Hence  V)W  varies  between  0-67  and  5.  Ratio  values  of  1,  3 and  9 were  selected  for  the  calculations. 

Let  T minutes  be  the  period  during  which  bearings  are  recorded  and  let  X nautical  miles 
be  the  distance  travelled  by  the  submarine  in  772  minutes.  The  values  of  P\ X are  recorded  below. 

(i)  Submarine  speed  = 5 knots 


(ii)  Submarine 


i 


Ratio  values  of  I,  7 and  49  were  selected  for  P/X. 

Using  selected  values  of  P and  V,  the  approximate  ATB  from  the  “bearings  only”  formula 
is  compared  to  the  known  value  of  the  ATB  for  a large  number  of  values  of  Ao  and  B.  E2  is 
the  “bearings  only"  estimate  of  the  ATB.  The  results  of  this  calculation  are  tabulated 
in  Appendix  II. 

In  the  tabulation  of  results,  Z represents  the  absolute  error  which  would  be  found  if  the 
bearing  lines  corresponding  to  A0  and  £2  were  drawn.  All  such  errors  > 20°  are  indicated  by 
an  asterisk  in  the  print  out  of  results.  An  examination  of  the  results  indicates  that  the  “bearings 
only”  estimate  for  ATB  improves  as  the  ratio  V/W  increases.  This  finding  is  also  reflected  in 
the  decreasing  incidence  of  asterisks  in  the  print  out  as  V/IV  increases.  The  error,  Z,  appears 
to  be  independent  of  the  ratio  P/W. 

CONCLUSION 

The  calculations  indicate  that  the  “bearings  only”  approximation  for  angle  to  bow  (ATB) 
on  the  convoy  is  critically  dependent  on  the  assumption  that  the  motion  of  the  submarine 

. . . , , „ , speed  of  convoy 

may  be  neglected.  The  greater  the  value  of  the  ratio ; — - — : — , the  smaller  the  error 

speed  of  submarine 

in  the  estimated  ATB  obtained  from  the  “bearings  only”  approximation.  The  error  in  the 
estimated  ATB  appears  to  be  independent  of  the  separation  between  the  submarine  and  the 
convoy. 
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APPENDIX  I 


DERIVATION  OF  APPROXIMATE  ANGLE  TO  BOW 

Case  ( i ) 

Application  of  the  Sine  Rule  to  triangles  OL'S,  OST  (See  figure  3)  yields  the  equations 

OS  US 


sin  (180°  - E2) 
OS 


sin  ( B — C) 
ST 


sin  0j  sin  (C  — E) 

Since  US  = ST  equation  (1)  may  be  simplified  as 

OS  ST 


sin  E2  sin  (5  — C) 

E2  is  an  exterior  angle  to  the  triangle  UOT.  Hence 

E2=01  + (B-C)  + (C- E) 
= e1  + b - £ 

Therefore  0t  = E2  — B + E 
Substituting  for  0X  in  equation  (2)  yields 

OS 


(1) 

(2) 

(3) 


ST 


sin  (£2  — B + E)  sin  (C  — E) 
Elimination  of  OS  and  ST  from  equations  (3)  and  (4)  leads  to 

sin  (E2  — B -j-  E)  sin  ( C — E) 
sin  E2  sin  ( B — C ) 

sin  EjCos  (B  — E)  — cos  E2  sin  (B  — E)  sin(C  — E) 


(4) 


sin  E2 

cos  (B  — E)  — 


tan  Ei  - 


cos  (fl  — E)  — 
sin  (B 


sin  ( B — C) 
sin  (B  — E)  sin(C  — E) 
tan  E2  sin  ( B — C) 
sin  (C  — E)  sin  ( B — E) 


(5) 


sin  (B  — C) 
C)  sin  (B  - E) 


tan  E » 


E2  = artan 


sin  ( B — C)  cos  (B  — E)  — sin  ( C — E) 
sin  ( B — C)  sin  (B  — E) 


(sir 
sin  (B  - C 


(6) 


C)  cos  (fl  ~ E)  — sin  (C 

Case  (//) 

Application  of  the  Sine  Rule  to  triangles  OLM,  OMN  (see  figure  5)  yields  the  equations 

OM  LM 


sin  02 
OM 


sin  (E  — C) 
MN 


sin  (E2  — 180°)  sin  (C  — B) 
Since  LM  — MN  equation  (8)  can  be  simplified  as 

OM  LM 

— sin  Et  sin  (C  — B) 


(7) 

(8) 

(9) 


8 


Since  360°  — Et  is  an  exterior  angle  to  triangle  OLN, 

360°  _ £2  = 0,  + (£  - C)  + (C  - B) 
= 0,  + E - B 


Therefore 

Hence 


02  = 360°  - £,  + B - E 
sin  08  = sin  (360°  — (Et  — B + £)) 
= — sin  (E2  — B + E) 


Equation  (7)  then  simplifies  to 

OA/ (10) 

- sin  (E2  — B + E)~  sin  (£  - C) 

Elimination  of  OM  and  LM  from  equations  (9)  and  (10)  leads  to 

sin  (£2  — B + E)  sin  (£-C) 
sin  E2  s>n  (C  — 

sin  £2  cos  (£  — £)  — cos  Et  sin  (B  - £)  _ sin  (C  - £) 

~ sin  £2  " ~ ~ sin  1 8 - c> 

This  equation  is  identical  with  equation  (5)  and  therefore  leads  to  the  formula  given  by 
equation  (6). 


APPENDIX  H 


THE  CALCULATOR  PROGRAM 

The  calculator  program  is  written  in  BASIC  for  the  Wang  2200A  Programmable  Calculator. 
A program  listing  is  provided  in  the  appendix.  Statements  180,  200  and  210  are  designed  to 
avoid  unnecessary  calculations  in  situations  where  the  bearings  only  approximation  is  not 
appropriate.  Statement  360  terminates  calculations  for  a given  value  of  B when  a collision 
course  is  indicated  by  the  artan  subroutine.  Statements  380,  390  and  400  terminate  calculations 
for  a given  value  of  B when  the  differences  between  B,  C and  E become  too  small  to  be  significant. 
Statements  420  to  550  terminate  calculations  for  a given  value  of  B when  calculations  indicate 
that  a collision  course  has  been  followed  by  ship  and  submarine.  Statements  560  to  610  check 
that  B^C^E  or  B<C<E.  If  neither  condition  is  met  the  calculation  is  terminated  for 

the  particular  value  of  B.  . 

The  artan  subroutine  is  contained  in  statements  790  to  990.  Limited  printout  is  produced 
when  “collision”  conditions  are  indicated.  A more  complete  data  printout  is  provided  for  each 
case  in  which  computations  are  completed.  The  subroutine  in  lines  1000  to  1040  controls  the 
number  of  lines  of  printing  per  page  of  output  and  arranges  page  numbering. 

In  the  printout  of  results  all  errors  > 20°  are  indicated  by*. 


:| L_  J| 

rr 


i 


10  REM  "BOR  1"  F7  TAPE  GJM4  18.  8.  1976 
20  SELECT  PRINT  215 <805 
30  SELECT  D 
40  L5=l 

50  PRINT  : PRINT  .PRINTUSING  1050,  L5 : PRINT  : PRINT 

60  PRINTUSING  1060: PRINT 

70  Ll=50 

80  L2=0.  000001 

90  L3=0.  000001 

100  L=0:fl*=,,*,‘ 

110  W=1 

120  FOR  I=1T0  3 

130  V=3~< I-1>*W 

140  FOR  J=1T0  3 

150  P=7^ < J-l ) ♦W 

160  FOR  R0=0  TO  340  STEP  20 

170  FOR  B=0  TO  180  STEP  20 

180  IF  (A0+BK10  THEN  730 

190  IF  V>2  THEN  220 

200  IF  RBS < B— R0 > <L2  THEN  730 

210  IF  RBS < B+R0— 180 ) <L2  THEN  730 

220  X1=V*SIN<A0> 

230  Y1=P+V*COS<R0> 

240  X2=2*X1 

250  V2=P+2*V*COS<A0> 

260  X3=W*SIN<B) 

270  V3=W*C0S<B> 

280  X4=2*X3 

29©  V4=2*V3 

300  X9=X1-X3:  Y9=Y1-Y3 

310  19=1 

320  GOSUB  790 

330  C=El+B-90 

340  X9=X2-X4:  V9=V2-V4 

350  GOSUB  790 

360  IF  I9<0  THEN  730 

370  E=El+B-90 

380  IF  ABS<B-EXL3  THEN  730 

390  IF  RBS < B > +RBS < C ) +RBS < E > <L3  THEN  410 

400  IF  RBS < B-180 > +RBS < C— 180 > +RBS ( E— 180 ) >L3  THEN  420 

410  PRINTUSING  1070,  V,  P,  A0,  B,  C,  E,  B-C,  C-E  . GOSUB  1000:  GOTO  730 

420  IF  RBSCBXL3  THEN  510 

430  IF  RBS < B-180 XL3  THEN  540 

440  IF  RBS<CXL3  THEN  480 

450  IF  RBS < C-180 ) >L3  THEN  560 

460  IF  RBS<EXL3  THEN  490 

470  GOTO  560 

480  IF  ABS<E-180>>L3  THEN  560 

490  PRINTUSING  1070,  V,  P,  A0,  B,  C,  E,  B-C,  C-E  : GOSUB  1000 
500  PRINT  TAB <7>i  "COLLISION" : GOSUB  1000:  GOTO  730 
510  IF  RBS < C-180 XL3  THEN  490 
520  IF  ABSCE-180XL3  THEN  490 
530  GOTO  560 

540  IF  flBS<CXL3  THEN  490 


J 


550  IF  flBS<EXL3  THEN  490 
560  10=0 

570  IF  B— C<0  THEN  580:  10=10+1:  GOTO  590 
580  10=10-1 

590  IF  C-EC0  THEN  600:  10=10+1:  GOTO  610 
600  10=10-1 

610  IF  RBSCI0X1  THEN  730 

620  E0=S I N < B-C  > *COS  < B-E  > -S I N < C-E  > 

630  IF  RBS<E0>>=L3  THEN  640 : E2=90 : GOTO  670 
640  E9=SIN<B-C>*SIN<B-E>/E0 
650  E2=RRCTRN<E9> 

660  IF  E2>0  THEN  670 : E2=E2+180 

670  IF  I0>1  THEN  680:  E2=E2+180 

680  2=RBS<E2-R0> : IF  2<180  THEN  690:  2=2-;fc«0 

690  IF  BBS <2 ><20  THEN  710:  PRINTUSING  1070,  V,  P,  R0,  B,  C,  E,  B-C,  C-E,  E2,  2,  R* 

700  GOTO  720 

710  PRINTUSING  1070,  V,  P,  R0,  B,  C,  E,  B-C,  C-E,  E2,  2 

720  GOSUB  1000 

730  NEXT  B 

740  NEXT  R0 

750  NEXT  J 

760  NEXT  I 

770  PRINT  HEXC 07 > 

780  END 

790  IF  RBSCX9XL3  THEN  930 
800  IF  RBSCV9XL3  THEN  900 
810  R=SQR<X9"2+V9"2> 

820  S0=SGN  < V9/R  > 

830  T1=V9/X9 
840  E1=RRCTRN  < T1 > 

850  IF  T1>0  THEN  880 

860  IF  S0<0  THEN  980 

870  E1=E1+180 . GOTO  980 

880  IF  S0>0  THEN  980 

890  E1=E1-180 : GOTO  980 

900  IF  X9>0  THEN  920 

910  El=180 : GOTO  980 

920  E1=0 : GOTO  980 

930  IF  RBS< V9>>=L3  THEN  950 

940  PRINT  "COLLISION”:  I9=-l:  GOSUB  1000:  RETURN 
950  IF  V9>0  THEN  970 
960  El=-90 : GOTO  980 
970  El=90 

980  IF  E1>0  THEN  990:  El=El+360 
990  RETURN 
1000  L=L+1 

1010  IF  L<L1  THEN  1040:  PRINT  HEX<0C> : PRINT  : PRINT 
1020  L5=L5+1 : PRINTUSING  1050,  L5 : PRINT  : PRINT 
1030  PRINTUSING  1060:  PRINT  : L=0 
1040  RETURN 

1050X  ### 

1060X  V P R0  B C E B-C  C-E  E2  2 

1070X  # ###  ###  ###  -###.  ##  -###.  ##  -###.  ##  -###.  ##  -###.  ##  -###.  ##  # 
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